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In this paper we present a digital holographic microscope to track gold colloids 
in three dimensions. We report observations of 100 nm gold particles in motion in 
water. The expected signal and the chosen method of reconstruction are described. 
We also discuss about how to implement the numerical calculation to reach 
real-time 3D tracking. @ 2012 Optical Society of America 



1. Introduction 

The tracking of nano-scale biological markers in motion is a challenging project, that could lead 
to a better understanding of dynamics in biological systems. Since digital holography gives access 
to the whole light field in a volume from a single hologram [1,2], this technique is suitable to 
reconstruct the field scattered by gold particles in three dimensions, and then determine their 
positions. This localization can be done without any mecanical scanning of the setup. The 3D 
tracking in digital holography was demonstated for micron-size particles in several publications 
with in-line holography [3-6]. The observation of nano-scale gold particles was also demonstrated 
in off-axis holography by Atlan et al. [7], and on the membrane of a living cell by Joud et al. [8]. 
Nano particles could be observed with the combination of off-axis geometry and phase-shifting 
technique [9], that achieves almost the theoretical noise level [10]. In this paper, we present first 
our digital holographic microscope, then we explain the reconstruction method. Some results of 3D 
tracking of 100 nm nano-beads are presented. 
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Fig. 1. Experimental setup. 



2. Experimental setup 

We chose to track lOOnm gold beads in Brownian motion in water. This particles are good light 
scatterers because their optical index in visible light is very low. Moreover they are non cyto- 
toxic and chemically inert [11]. The sample is prepared using a 400 /jm deep chamber between 
two cover glasses filled with a nano-beads solution. The experimental setup is shown on the figure 
1. The optical setup use a 100 mW DPSS laser emitting at 660 nm. The laser beam is splitted 
with a polarizing beam splitter. We adjust the ratio of energy between both beams with an half- 
wave plate placed before the splitter. The first beam, called illumination beam, is focalised on the 
sample. The second beam, called reference beam, is spatially filtered and expanded in a way to 
cover all the CCD chip. A dark-field microscope objective from Edmund Optic (NA = 0.5 and 
36 x magnification) collects the scattered light and blocks the illumination beam. This dark field 
configuration is necessary to avoid the saturation of the CCD chip : the Rayleigh-Mie scattering 
model [12] gives a cross section of 0.015 /xm 2 for a 100 nm diameter gold particle at 660 nm. 
During 1 ms, the number of incident photons is typically 1014 photons at 660 nm, and the area 
of the illumination beams waist is ~ 4 x 10 -8 m 2 . The number of scattered photons by the nano- 
bead collected through a numerical aperture of 0.5 is equal to 8 x 10 6 photon during 1 ms. The 
ratio of number of photons between the illumination beam and the scattered light is about 10 8 , 
so its not possible to record the scattered signal on our 16 bits CCD chip without blocking the 
illumination beam. The dark field objective allows to use the full dynamic of the CCD chip to 
record the low scattered signal. A non-polarizing beam splitter behind the microscope objective 
combines the reference and the scattered signal. This splitter is titled of a few degrees to be in 
off-axis geometry [13]. The interference pattern is recorded on the 16 bits CCD chip of 1024 x 1024 
pixels. 
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Fig. 2. Steps of the numerical reconstruction. 



The observation of particles in motion imposes to limit the exposition time in order to limit 
the blurring of the signal. The magnification calibration of our setup gives a lateral pixel size of 
160 nm. In our case, the diffusion constant of the Brownian motion is equal to D ~ 4/im 2 .s _1 
according to the Stokes-Enstein law. So the exposition time must be lower than 3 ms to avoid that 
the signal is shiffted by one pixel. The exposition time of our camera can be set under this limit, 
but the acquisition period during 2 frames is limited to 80 ms. The travel of 100 nm particles in 
water during this period is too long to use the phase-shifting technique in order to improve the 
signal-to- noise ratio [9,14]. 

3. 3D reconstruction of the scattered light 

The camera record for each frame an intensity I cc d described in the following equation : 

Iced = Iref + L + ErefE* s + E s E* e j (1) 

where I re f and I s are respectively the intensity of the reference and the scattering light. E re t and 
E s are the respective fields. The tilted beam splitter which combines the two fields adds a spatial 
frequency on the cross terms of interference (the third and fourth terms of the equation (1)). The 
effect is that these two terms can be separated in the Fourier space of the holograms. Figure 2 
summarizes the steps of our reconstruction. First we multiply the hologram by a complex phase 
matrix exp[j-7r(x2 + y2)/{\d)\ where A is the wavelength and d the distance between the CCD and 
the output pupil of the microscope objective. This operation is similar to a back-propagation of 
the hologram to the output plane of the objective. This propagation compensates the sphericity 
of the wave added by the objective. Then we calculate the Fast-Fourier Transform (FFT) of the 
corrected hologram. In the Fourier plane, the two first term of equation 1 are in the center, while 
the cross terms are respectively centered on the frequency added by the off-axis geometry and its 
complex conjugate. We numerically filter the term E s E* e j of equation 1 with a circular mask of 
radius equal to the apparent radius of the output pupil. The filtered map of the spatial frequencies 
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(f[k x , ky]) obtained is multiplied by the following matrix of propagation : 

( i n x Sz).jX(kl +kl)\ 
K[k x , k y , n\ = exp I — y — I (2) 

where bz is the arbitrary step of axial propagation (typically 100 nm) and n an integer. Then we 
calculate the inverse FFT. The result e(n) = FFT -1 [/if] is the reconstruction of the field at the 
distance n x bz. We repeat this multiplication by varying the integer n to reconstruct the volume. 
We get a stack of lateral sections e(nbz) spaced by bz. The positions of the scatterers are determined 
by the local maximum of the intensity of the whole reconstructed field. 

All this numerical reconstruction is implemented on a graphic card unit. We developp our algoritm 
with the Nvidia CUDA architecture to parallelize all the calculation on the 448 cores of a Nvidia 
Geforce GTX 470 card. All the steps of the reconstructions, except for the FFTs, need multiplication 
pixels to pixels and not matrix multiplication. This kind of operation are perfectly adapted for 
parallelization on graphic cores [15-18]. The FFT steps are made with the CUFFT library. We 
divided the time of reconstruction by a factor 30 compare to the same reconstruction made only 
on the CPU (double quad-core Intel Xeon E5520 at 2.27 GHz). We developped a software for live 
visualisation of reconstruction using OpenGL library to display the datas calculated on the graphic 
card. 

4. Current results 

We are able now to observe lOOnm gold nano particles in Brownian motion during few tens of 
positions. After that we generally lose the targets. The figure 3 shows the observation we are able 
to do in the present state of the setup. We see on the figure the lateral section where the intensity 
of the scattered light by the nano-bead is maximum. The PSF of the reconstructed field is plotted 
for X, Y and Z directions. The lateral PSF width is equal to 1 /jm, and the axial width is equal to 
14 /mi. The signal-to-noise ratio of the reconstruction is better than 100. These results are obtained 
with an exposition time of only 1ms and an acquisition period of 80 ms between two frames. 

5. Conclusion 

We report in this paper the present performance of our digital holographic microscope. It consists 
of a Mach-Zender setup in off-axis geometry. A dark field microscope objective collects the light 
scattered by nano-sized gold particles, which interferes with the reference plane wave on the CCD 
camera. The scattered light field is numerically propagated from the interference pattern. We used 
our own algorithm based on the CUDA architecture to reduce the calculation time. The improve- 
ment clearly shows that this technology is a very promising tool for digital holography. We are able 
now to localize 100 nm gold nano particles in Brownian motion from holograms exposed during 1 
ms. We expect to obtain similar performance for gold particles in a crowed medium or in a living 
cell in order to study the diffusion law in these specific media. 
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Fig. 3. Intensity of the scattered field by a 100 nm gold particle plotted in the section where 
the intensity is maximum. The three curves at the left shows the PSF respectively in X, Y and Z 
directions. 
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